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Abstract 
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I. INTRODUCTION 



There has been a continues interest on the solutions of the Klein-Gordon (KG) equation 
and Dirac equation for some certain potentials. The KG equation is solved by assuming 
that the scalar potential equals or not equals to the vector potential, such as Rosen-Morse- 
type potentials, Morse- like potential [1-8]. On the other hand, it is also solved for the case 
of mixing of scalar and vector potentials, such as vector-scalar Coulomb, kink-like poten- 
tials, harmonic oscillator, Hartman potential, and Hulthen-type potential [9-13]. Further, 
the Dirac equation is solved for different types of potentials, such as harmonic, and linear 
potentials [14], for an uniform magnetic field [15], generalized asymmetrical Hartmann po- 
tentials [16]. Various methods are used in the solutions, such as Nikiforov-Uvarov (NU) 
method, by using the hyper geometric type equations, and separation of variables [17-21]. 

In recent years, the effects of the position-dependent the mass on the energy spectra 
and corresponding eigenvalues of the above equations has been received a great attention 
[22, 23]. The exact solutions of the above equations, and the Schrodinger equation in the 
context of spatially dependent mass have been studied by many authors [24-28]. Many 
different types of mass distributions have been used in literature, such as polynomial mass 
function, exponential, and hyperbolic mass distributions [30]. Our mass function is similar 
to that used in Ref [24] by Dutra, and Almeida. In the first example, the mass function, 
and potantial are taken as In the present work, we prefer an exponential type mass function 
to fnd out the energy spectrum and corresponding wave functions of the Hulthen potential. 
The Hulthen potential has applications in the wide range of different areas such as nuclear, 
and particle physics, atomic physics, solid state, and chemical physics [7]. 

In this study, we solve the radial KG equation for the Hulthen potential by using NU- 
method within the framework of an approximation to the centrifugal potential. The NU- 
method is based on solving the second order equation by reducing to a generalized equation 
of hypergeometric type [31]. 

The organization of this work is as follows. In Section II, we solve the radial part of 
the KG equation by using an approximation for the centrifugal term [33, 34], and compute 
energy eigenvalues of the bound states and the corresponding eigenf unctions. We give our 
conclusions in Section III. 
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II. NIKIFOROV-UVAROV METHOD AND BOUND STATES 



The radial part of the KG equation reads [29] 



(1) 



where ^ is the angular-momentum quantum number, E is the energy of the particle, c is the 
velocity of the light. 

The Hulthen potential has the form [32] 



= -V. . (2) 

where Vb, and (3 = l/a are constant parameters. 

Eq. (1) can not be solved exactly because of the centrifugal potential. So, one has to use 
an approximation for this term. This approximation can be taken as [33, 34] 



There are very different mass-distributions are used in the literature, such as an expo- 
nentially, and quadratically mass functions [24], trigonometric mass-distributions, and mass 
functions of the forms m(r) = r'', especially for three-dimensional problems [26]. Here, we 
prefer to use the following spatially dependent mass 



m(r) = mo - ^ _ , (mo > mi) . (4) 

which corresponds to a decay particle, and provides us an approximate solution of the radial 
part of the effective KG equation, mo and mi in this distribution are two arbitrary positive 
parameters. The mass function of that form enable us to check out the results in the 
limit of the constant mass. 2. It is well known that the Schrodinger equation (SE) should 
be reproduced in the case of position-dependent mass, because of the ordering- ambiguity 
problem between momentum, and mass operators in kinetic energy term. The kinetic energy 
term should be written as p 2m{x) P- When the mass depends on coordinate, and it can be 
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seen that the operators no longer commute[see for details Ref [10] by Gonul, and Ref [24]]. 
The ordering- ambiguity problem does not appear in the case of KG equation. 
Substituting Eqs. (2) and (4) into Eq. (1), and by using Eq. (3), we get 



-/3r 



-I3r]2 



+ 



2^2 



h'6 



-I3r 



+ 2EVo 



-l3r 



2 4 2momic 

-TTIqC + 



2 4 



1 - e-f^^ (1 - e-/3'-)2 
By using the transformation z — 1 — e~^'^ (0 < 2; < 1), we obtain 



1 - e-/3-)2 
0(r) = . 



(5) 



^ d(j){z) 1 r 2 2 2 21 , / N n fR\ 



where 



-al ^ Q\E'' - 2EVo + - mlc") , 

-al = Q2(2momiC^ + 2EVq - 2V^) + l{l + 1) , 

= Q\V^ - mY) - t{t + 1) . (7) 

and = 1/ pi^U^c^. Now to apply the NU-method [31], we rewrite Eq. (6) in the following 
form 



r(.-)+Ii|^'(.-)+^^W^O, (8) 

where (7{z) and d{z) are polynomials with second-degree, at most, and f{z) is a polynomial 
with first-degree. We define a transformation for the total wave function as 



<t>{z) = az)m- (9) 

Thus Eq. (8) is reduced to a hypergeometric type equation 



a{z)i)"{z) + T{z)i)'{z) + \i){z) = 0. 
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(10) 



We also define the new eigenvalue for the Eq. (8) as 



A = Xn^-nr'- !^i!i-ila",(n = 0,l,2,...) (11) 



where 



r{z) = f{z) + 2TT{z). (12) 

The derivative of t{z) must be negative. A(A„) is obtained from a particular solution of the 
polynomial ipn{z) with the degree of n. ipni^) is the hypergeometric type function whose 
solutions are given by [31] 



where the weight function p{z) satisfies the equation 



^Ja{z)p{z)]^T{z)p{z). (14) 
On the other hand, the function ^{z) satisfies the relation 



e{z)/az)=7r{z)/a{z). (15) 
Comparing Eq. (6) with Eq. (8), we have 

f{z) — —z , a{z) = z{l — z) , a{z) = —alz"^ — a\z — a\ (16) 
The t:{z) has the form 



, , a'iz)-fiz) , .a'(z) - f(z).„ ~/ v , / v 



or, explicitly 
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TT U 



(18) 



The constant k is determined by imposing a condition such that the discriminant under 
the square root should be zero. The roots of k are ki^2 = — ~ 2a| =F + 4a|, where 



A = ^/al + al + af. Substituting these values into Eq.(18), we get for 7r{z) 



n{z){k^k,)= ^{l-z)T 



and 



7r{z){k 



A- 



z + 



k,)^l{l-z)T[{A + ^[Y+^l)z-^ 



(19) 



(20) 



Now we calculate the polynomial r{z) from 7r(z) such that its derivative with respect to z 
must be negative. Thus we take the first choice 



T(z) = l-2^i +ai -2( 
The constant X — k + 7r'{z) becomes 



A 



- +ai +l]z. 



(21) 



and Eq. (11) gives us 



2al + A^l + 4al-l-[A-^\+al), 



(22) 



A„ = 2n(^A- +a| + ij -n. 



(23) 



Substituting the values of the parameters given by Eq. (7), and setting A = A„, one can find 
the energy eigenvalues as 



Enf — 



Vo 



+ 



SQ^miC^lrrii — 2mo) 



+ 8g^A^^(2moC^ - 2momic^ + m\c^ - V^) - ^'^^ 



(24) 
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where 



(2n+l) + -^l + 4ai . (25) 

We see that the energy levels for particles and antiparticles are symmetric, and the ground 
state energy is different from zero. 

We also get the energy eigenfunctions for the constant mass case for s-states 

K,,=o - Y ± ^ y AQHfTN^ ~ 160^' ^^^^ 

where 



A^' = (2n + 1) + 1 - AQ^Vo^ . (27) 

It is the same result with those in literature for ^ = 1 [21]. 

Now let us find the eigenfunctions. We first compute the weight function from Eqs. (12) 
and (14) 



p{z) = (1 - zf^ , (28) 

and the wave functions become 



^^l+4a2 ^-^ _ ^-j2A dz' 

where 6„ is a normalization constant. The polynomial solutions can be written in terms of 
the Jacobi polynomials [35] 



M^)=bnPi''^"'^\l-2z), 2A>-1, ^l + 4ai > -1 . (30) 
On the other hand, the other part of the wave function is obtained from the Eq. (15) as 



^n+y/ l+4a§ j-^ _ 2;)"'+2^ 



(29) 



(31) 



Thus, the total eigenfunctions take 



Uz) -b'Al- z ^ i'+V^] p(^^'^^) (1 _ 2z) . (32) 

where b'^ is the new normahzation constant. 

Total wave function for the constant mass case for any ^-state takes 



It is consistent with the results obtained in the literature [21]. The parameters for this case 
are given by A' — •\/a^~-i-aJ~+a^, and 



4^ = 2Q^VoiVo-E) -£{£+!), 

a'i = -QVo' + £(£+1). (34) 

III. CONCLUSION 

We have obtained approximate solution of the radial part of the effective mass KG equa- 
tion for the Hulthen potential in the framework of an approximation to the centrifugal term 
for any £ values. We have obtained the energy spectra and the corresponding radial part 
of the wave functions by applying the NU-mcthod. Wc have found a real energy spectra 
for the Hulthen potential in the PDM case. To check our results, we have also calculated 
the energy eigenvalues and eigenfunctions of the particle and antiparticles for the case of 
constant mass limit for s-waves, and seen that the results are consistent with those in the 
literature [21]. 
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